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Abstract 

We previously showed that the real eigenvalues of 3 x 3 octonionic Hermitian ma- 
trices form two separate families, each containing 3 eigenvalues, and each leading to an 
orthonormal decomposition of the identity matrix, which would normally correspond 
to an orthonormal basis. We show here that it nevertheless takes both families in order 
to decompose an arbitrary vector into components, each of which is an eigenvector of 
the original matrix; each vector therefore has 6 components, rather than 3. 



1 Introduction 

In previous work [p], @, [J we considered the real eigenvalue problem for 3x3 octonionic Her- 
mitian matrices. [] For this case, there are 6, rather than 3, real eigenvalues ||. We showed 
that these come in 2 independent families, each consisting of 3 real eigenvalues, which satisfy 
a modified characteristic equation rather than the usual one. Furthermore, the correspond- 
ing eigenvectors are not orthogonal in the usual sense, but do satisfy a generalized notion of 

*Supported in part by U. S. Department of Energy Grant No. DE-F402-ER 40685. 

1 There are also eigenvalues which are not real [Q, g|, whose intriguing physical properties were discussed 
in |g, Q]. A related eigenvalue problem, which admits only real eigenvalues, was discussed in ||. 
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orthogonality. Finally, all such matrices admit a decomposition in terms of (the "squares" 
of) orthonormal eigenvectors, although these matrices are not idempotents (matrices which 
square to themselves). 

It is the purpose of this paper to extend these results by showing how to project an 
arbitrary 3-component octonionic vector into the eigenspaces determined by the real eigen- 
values of a given 3x3 octonionic Hermitian matrix. Somewhat surprisingly, this turns out 
to require all 6 eigenvalues; generic 3-component octonionic vectors therefore have 6 vector 
components, rather than 3. 



2 Octonionic Eigenvectors 

In this section, we summarize the main results of our previous work |], ^ 

The (real) octonions O are the noncommutative, nonassociative division algebra gen- 
erated by adjoining 7 anticommuting square roots of — 1, with a particular multiplication 
table, to the real numbers M. Q This generalizes both the complex numbers C and Hamil- 
ton's quaternions H; further details can be found in [0 or [[H^, 11]. We refer to a matrix as 
complex or quaternionic if its components lie in some complex or quaternionic subalgebra 
of O, respectively. 

Let A be a 3 x 3 octonionic Hermitian matrix. We consider solutions of the eigenvalue 
equation 

Av = Xv (1) 
where v G O 3 and A Gl. Then A satisfies a modified characteristic equation of the form 

det(AJ - A) = A 3 - (tr A) A 2 + a{A) A - det A = r (2) 

where the determinant is defined in terms of the Freudenthal product and where 

a(A) = i((tr^) 2 -tr(A 2 )) (3) 

The real number r is a root of the quadratic equation 

r 2 + 4$r - H 2 = (4) 

where both the associative 3-form $ and the associator a are totally antisymmetric functions 
of the non-real components of A, which are given explicitly in Section |J There are 3 real 
solutions for A corresponding to each solution for r, which labels the families. 

If the roots r\ and r2 of (f|) are the same, then, over the real octonions, $ = = a, 
which forces A to be quaternionic (and ri = = r?). Unless A is complex, there are still 
two families of eigenvectors, but r can not be used to label them. Rather, as described in 
the next section, one family consists of the usual quaternionic eigenvectors, while the other 
eigenvectors are purely octonionic, in the sense that their components are orthogonal (as 



2 This construction can be performed over any field K, although this will not in general yield a division 
algebra. We restrict ourselves here to the real octonions, as was done in pi 0. As was explicitly verified 
in H , most of these results hold over an arbitrary field, since division is only used in constructing normalized 
eigenvectors. 
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vectors in R 8 ) to the quaternionic subalgebra containing the components of A. Finally, if 
A is complex, there is only one family of eigenvectors, obtained by (right) multiplying the 
usual complex eigenvectors by arbitrary octonions. 

Within each family, many of the usual properties hold. Two eigenvectors u, v corre- 
sponding to different eigenvalues in the same family are orthogonal, but in the generalized 
sense that 

(uw)v = = (vv*)u (5) 

Even in the case of repeated eigenvalues, 3 normalized, orthogonal eigenvectors v m , with 
eigenvalues A m , can be found in each family. Orthonormality can be expressed in the form 

3 

m=l 

where X denotes the identity matrix, and these eigenvectors can be used to decompose A in 
the form 

3 

-4=5^ ^ m v m vl (7) 

m=l 



3 Quaternionic Projections 



In this section, we consider the much easier case of quaternionic matrices. While we are 
primarily interested in the 3x3 case, all results in this section apply tonxn quaternionic 
Hermitian matrices. 

It is easily shown |] that the right eigenvalue problem 

Av = v\ (8) 

for quaternionic Hermitian matrices A admits solutions only for A £ M. It is straightforward 
to find an orthonormal basis {v n } of (quaternionic) eigenvectors, which satisfies 

Av m A m f m (9) 

V ln V n = 6 mn (10) 

Given any quaternionic vectors v, x, the projection of x along v is simply v(v'x), so that 
any quaternionic vector has vector components of the form 

x = ^2v m (vl n x) (11) 

which makes sense, since (the parentheses are not needed and) orthonormality can be written 

I= ^2v m vl (12) 



as 



We now ask about the octonionic properties of the quaternionic matrix A. To do this, 
we view the octonions through the Cayley-Dickson process as 

o = e©ffl (13) 
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where £ is any octonionic unit orthogonal (as vectors in R 8 ) to all elements of H, the quater- 
nionic subalgebra of O containing the components of A. 

What are the octonionic eigenvectors of A! As discussed in jj], if the components of A 
and v lie in the (same) quaternionic subalgebra H, then we have 

A{£v) = £{Av) (14) 

Thus, if u is an eigenvector of A, then lu is an eigenvector of A with the same eigenvalue! 
This implies that, over O, the real eigenvalues of A consist of its eigenvalues over HI together 
with the eigenvalues of A over HI. 

As expected, there are two families of (octonionic) eigenvectors of A, the usual family 
from the quaternionic eigenvalue problem, together with an "extra" family, namely {£u m } 
where {u m } are the (quaternionic) eigenvectors of A, with real eigenvalues [i m . This latter 
family also leads to a decomposition of A. To see this, note that if u is quaternionic then 

{£u){£ U y = £{uv))J = W (15) 

from which it follows that for any quaternionic u, x 

[(£u)(£ U y] (£x) = ^J{£x) = £(uu*x) (16) 

But this means that (£u)(£u)^ has precisely the same projection properties over £M as uv) 
has over HI. Thus, if {u m } are (quaternionic) orthonormal eigenvectors of A, then {£u m } are 
(octonionic) orthonormal eigenvectors of A, and in particular 

Z = X>«J(&4) (17) 
A = J2^u m )(£ul) (18) 

Putting this all together, we can decompose an arbitrary octonionic vector x in terms 
of the eigenvectors of A and A. First of all, there is a unique decomposition of x into a 
quaternionic piece and a "purely octonionic" piece, given by 

x = x\ + £x2 (19) 

where both x\ and X2 are quaternionic. The rest is easy: Expand x\ in terms of the (quater- 
nionic) eigenvectors of A, and Xi in terms of the (quaternionic) eigenvectors of A. Equiva- 
lently, expand x\ in terms of {v m } and £x<i in terms of {£u m }. This leads to 

X = ^2 v m( v L x l) + Yl {( iu m)( iu m)) ( ix ^ 

= Yl V m V l X l + 1 Yl U m U L X 2 (20) 

We have therefore succeeded in projecting x into vector components, each of which is an 
eigenvector of A. However, there are twice as many components as one would have expected 
from the purely quaternionic problem. Furthermore, while it is clear that v m ^xi is an 
eigenvector of A with eigenvalue A m , since it is a quaternionic right multiple of v m , it is 
rather surprising that f^UmU^Xi) is an eigenvector of A, since it is not in general a multiple 
of £u m \ 

We will see in the next section that this entire structure carries over to the fully octonionic 
setting, but that the identification of the two families is not as easy. 

3 If A is complex, the choice of HI C O is not fixed. In this case, however, there is only one family of 
(octonionic) eigenvectors. We assume in what follows that A is not complex. 
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4 Decompositions of the Octonions 



In this section, we summarize and interpret some of the more technical results from [|J, 
which we will need to prove our main theorems in Section |8|. These results are interesting 
in their own right, however, as they lead to the consideration of several intriguing ways of 
decomposing the octonions into subspaces. 

Write the components of the octonionic Hermitian matrix A as 




A=\a e c (21) 



with d,e, f G K and a, b, c G O. With these conventions, the associative 3-form 0] is given 
by 

$ = Re(a xbx c) = -Re(a(6c) - c(ba)) (22) 
and the associator is given by 

a = [a, b, c] = (ab)c — a(bc) (23) 

Viewing the octonions as the real vector space M 8 , we introduce the subspace TcO spanned 
by the elements of A. In other words, T is the real vector space consisting of all (real) linear 
combinations of the components of A, that is 

T=(l,a,b,c) (24) 

Since the quaternionic case was covered in the previous section, we will assume throughout 
the remainder of the paper that a ^ 0, which in turn implies that T has dimension 4 (as a 
real vector space). 

Define the characteristic operator K associated with A by 

K[x] = A{A{Ax)) - {trA){A{Ax)) + a(A){Ax) - (detA)x (25) 

where x G O 3 . As we may readily verify, the operator K is diagonal, and can thus be 
reinterpreted as an operator acting on single octonions. This operator, also called K, can be 
written in the form 

K[p] = c{b(ap)) + a(b(cp)) - (c(ba) + a{bc)^jp (26) 

for any p G O. In this form, K is the same as the operator introduced in ||, where it was 
shown to satisfy the same quadratic equation as r, that is (compare (Q)) 

K 2 + — \a\ 2 = (27) 

In fact, comparing ( [25] ) with (g), we see that if v is an eigenvector of A, so that ([!]) holds, 
then 

K[v] = rv (28) 
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with r a solution of (|J). Thus, the operator K can be used to distinguish the 2 families of 
eigenvectors, labeled by the 2 possible values of r. 

Using results from [[| , we can show that K takes a particularly simple form on T, namely 

K[t] =ta (te T) (29) 

Since 

T 1 - = Ta (30) 

it follows using (|27|) that 

K[u] = -u(a + 4$) (tie? 1 ) (31) 

An intriguing property of the decomposition 

O = T © Ta (32) 

is that it is almost of the Cayley-Dickson form. Direct computation establishes the multi- 
plication table (for ti,t 2 G T) 

ti(t 2 a) = (t 2 ti)a (33) 
{t x a)t 2 = {t x T 2 )a (34) 
(tia)(t 2 a) = — t 2 ti\a\ 2 (35) 

which is of the Cayley-Dickson form. However, this is not a true Cayley-Dickson extension 
of T, since T is not closed under multiplication. 

Using these ideas, it can be shown (Proposition 3.1 of |3|]) that the eigenspace of K with 
eigenvalue r is precisely T(r + 4$ + a), that is 

K[q] = r m q 3t G T : q = t(r m + 4$ + a) (36) 

where the solutions of (|]) are denoted r m , with m — 1, 2. If we define Q 

r m + 4$ + a 

Sm ~ 2(r m + 2$) {67) 
then the decomposition of O into its eigenspaces under K takes the form 

O = Tsi © Ts 2 (38) 
and the normalization is chosen such that 

si + s 2 = 1 (39) 



4 Recall that, over the real octonions, the two solutions of (Q) are distinct unless $ = = a, in which 
case A is quaternionic and can be handled by the methods of the previous section, and K is not needed to 
distinguish the two families. 

J Since r\ + + 4$ = 0, the denominator can only vanish if r\ = r%, or equivalently r = in the notation 
of ||; this cannot happen here since we are assuming a^O. 
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The eigenspaces T m := Ts m are orthogonal, since K is self-adjoint. Explicitly, the inner 
product in IR 8 takes the form 

2p ■ q = pq + qp = pq + qp (p,q G O) (40) 

and K satisfies 

K[p]-q = p-K[q] (41) 
It is now straightforward to define projection operators 

= K + r m + A<S> 
Km 2(r m + 2$) 

from O to T m , which satisfy 

K x + K 2 = 1 KxK 2 = = K 2 K 1 K\ = K x K\ = K 2 (43) 

Finally, since ~s~[s 2 is proportional to a, we also have 

T 2 = Tia (44) 

The following two results are key to understanding the family structure. First, assertions 
(v) and (vi) of Proposition 3.2 of imply 

p,qeT m =}► WGT (45) 

Second, as we also see from Proposition 4.1 of ||, for any p±,p 2 G T and g G T m , we have 

\pi,P2,q\ =Pzq (46) 

where p 3 G O depends on p x , p 2 (and m) , but does noi depend on q. These two results are 
the primary ingredients in the proofs in Section p] of our main results. 



5 Octonionic Projections 

We now return to the discussion of octonionic eigenvectors, showing how to project any 
vector into the 6 eigenspaces they determine. We begin with the following result. 

Theorem 1: Let a 3 x 3 octonionic Hermitian matrix A be given, and suppose that 
v G O 3 is an eigenvector of A, that is, Av = Xv with A el. Suppose further that K[v] = rv, 
with K as in fi2f\), and let y G O 3 be any other vector satisfying K[y] = ry (for the same r). 
Then 

(vv r ) ((vv^y) = (v^v) ((vv^y) (47) 

Proof: See Section |^. 

In other words, (vv')y is an eigenvector of a new 3x3 octonionic Hermitian matrix, 
namely vv'\ if v is normalized, the eigenvalue is 1. Note, however, that (vv^)y is not neces- 
sarily a multiple of v\ If u is another eigenvector of A in the same family, so that K[u] = ru 
and Au = pu and where we assume A ^ jj G K, then u and v are orthogonal (in the 
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sense (|5|)). This in turn means that u and (vv^)y are eigenvectors of the matrix vv^ with 
different eigenvalues, and hence must themselves be orthogonal, that is 

uv) ((W)y) = (48) 

This shows that, restricted to the appropriate K eigenspace, uv) and vv* are projection oper- 
ators! Extending this to a family {u,v, w} of orthonormal eigenvectors, we obtain using (0) 

A ((vv*)y) = A ((vv*)y) (49) 

Furthermore, (§) leads to 

y = (uu*)y + (vv*)y + (ww')y (50) 

which expresses y as the sum of its projections along the eigenvectors. 

Given an arbitrary vector i 6 O 3 , we can use the projection operators K m to project x 
into the eigenspaces of K, by defining 

x n = K n [x\ (51) 

and noting that 

x = xi + x 2 (52) 

We can then use (0) to expand each of x m along the appropriate eigenvectors. Putting this 
all together, any vector x can be expanded into six terms, consisting first of the projections 
into the 2 eigenspaces of K, then decomposing each of these separately into 3 eigenspaces 
of A 



6 Families of Eigenvectors 

Is the family structure a property of the matrix A, or of its eigenvectors? We show in this 
section that it is the latter. 

Decompose A in terms of an orthonormal eigenbasis as in (|7|), yielding 

A = Xuu' + fj,vv' + vwvj^ (53) 

Then, by the results of the last section, all of u, v, w have the same eigenvalue under K, 
namely 

r = Xfiu - det A (54) 

Note that this result holds regardless of the values of A, fi, v\ it is a property of the orthonor- 
mal "basis" {u,v,w}, not of a particular matrix A. 

Consider the special case where A = vv', with v'v = 1. Then v is an eigenvector of A 
with eigenvalue 1, and 

K[v} = -(detA)v (55) 

since 

tvA=l a(A) = (56) 

In this case, fl4"T|) simplifies to 

K[y] = — (detvv')y <^=^ vv^ ((vv')yj = (vv')y (57) 
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since the converse follows immediately from the definition of K (with A = vv^). Special 
cases of vectors u satisfying this condition are u orthogonal to v, so that 

Au = (vv ] )u = (58) 

or 

uu 1 = svv^ 1 (59) 

with s G R, in which case 

Au = (vv^)u = -{uv))u = -u{v)u) = u (60) 

The decomposition ( |50|) shows that any y satisfying 

K[y] = -(detA)y (61) 

can be written as the sum of terms of one of these two types. 

If vv' is complex, then v determines a single "family" containing all vectors. Furthermore, 
v must then be a (right) octonionic multiple of a complex vector t> ; the "family" structure 
in this case corresponds to completing vq as usual to a complex orthonormal basis. However, 
there are also other families containing v. (If v is a real eigenvector of A, then A must be 
real or complex, and there is then no nontrivial family containing v.) 

If vv' is not complex, there is an 8-parameter family of vectors orthogonal to v, and a 
4-parameter family of vectors satisfying fl5"9]). We therefore see that v alone determines a 
12-parameter family of vectors satisfying both fl47Q and fl6T|) . This family consists precisely 
of the (real) linear combinations of a family of eigenvectors of any matrix of the form fl53p. 
In other words, a family of eigenvectors is just an orthonormal "basis", but each such basis 
spans only half of O 3 . 

We can complete this picture by showing that each vector v (with vv' not complex) lies in 
exactly one such family. To do this, we first remove the requirement that v be an eigenvector 
from (pT7|). Specifically, we have 

Theorem 2: Let K be the operator constructed via ([29i) from a given 3x3 octonionic 
Hermitian matrix A, and let y, z € O 3 . Suppose K[y] = ry and K[z] = rz with r£l. Then 

(yyl) ((yy^z) = (y^y) ((yy^)z) (62) 

Proof: See Section 

This finally allows us to argue that if y is in the family defined by v, by virtue of satisfying 
either (and hence both) of ([17]) and (|B~TD, then v is also in the family defined by y. This 
follows immediately from Theorem 2, since ( pT| ) is enough to satisfy the condition of the 
theorem, and (j62|) (with z = v) ensures that v is in the family defined by y. 

More generally, Theorem 2 shows that u, w belong to the same family (with u normalized 
by u^u = 1) if and only if 

uu* ((uu')w) = ((uu')w) (63) 

For, (|63| ) shows that w is in the family defined by u. If u is in the family defined by v , 
say, then the previous argument shows that v is also in the family defined by u. But then 
v and w are both in the family defined by u, and the theorem shows that they must each 
be in the family defined by the other. This shows that any family which contains u must 
also contain w. The converse, namely that u, w in the same family implies (|63|), follows 
immediately from the theorem. 
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7 Discussion 



As made clear in Section |], the family structure, which was originally discovered on O 3 
arising from the eigenvectors of A, can be viewed as a property of the octonions themselves. 
From this point of view, the only purpose of A is in determining the octonions a, b, c used 
to define T. 

So suppose a, b, c G O are given, and further assume that a = [a, b, c] ^ 0. Does this 
decomposition of O depend on the choice of a, b, c, or merely on T? Suppose A, B,C G T, 
so that 

A = A + A x a + A 2 b + A 3 c (64) 
with A m G 1R, and similarly for B and C. If we write A for (Ai, A 2 , A3) G M 3 , etc., then 

[A, B, C] = (A, B, C) [a, b, c] (65) 
$(A,B,C) = (A,B,C)$(a,b,e) (66) 

where (A, B, C) = (Ax B) ■ C is the vector triple product, assumed to be nonzero. Since both 
a and $ scale with the same factor, so does r — which means that the s m are unchanged! 
The decomposition of O thus depends only on the subspace T. 

In terms of matrices, this means that any matrix with non-real components A, B, C G T 
satisfying (^4, B, C) 7^ determines the same family structure as the original matrix A. The 
2 families determined by A are just 

O 3 = F!©F 2 := (T!) 3 ©(T 2 ) 3 (67) 

An essential ingredient in this picture is the result (f45D, which ensures that, for any vector 
v G F m , the components of vv* are in T! 

We conclude by emphasizing that the subspaces T m have many of the properties of the 
quaternions: They are 4- dimensional, they lead to an almost Cayley-Dickson multiplication 
table for the octonions, and, our main result, they have just enough associativity to enable 
arbitrary vectors to be projected into eigenspaces. 
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8 Proofs of the Theorems 

We give here the proofs of our two theorems. This section can be skipped on first reading; 
all results have also been independently derived using Mathematics. 

Theorem 1: Let a 3 x 3 octonionic Hermitian matrix A be given, and suppose that 
v G O 3 is an eigenvector of A, that is, Av = Xv with A el. Suppose further that K[v] = rv, 
with K as in (21.), and let y G O 3 be any other vector satisfying K[y] = ry (for the same r). 
Then 

(vv*) ((vv*)y) = (v*v) ((vv*)y) 



Proof: By (|36]), v,y G (T m ) 3 , and by (f45"D the components of vv* are in T. Thus, 
introducing a natural matrix associator, (fOT) implies that 



(vv*) ((vv*)y) = ((vv*)(vv*)) y + [vv*, vv', y] = Vy (68) 
where V does not depend on y\ But alternativity implies that the vector associator satisfies M 

[v,v*,v}=0 (69) 

which in turn implies that 



(vv 



^ ((vv*)v) = (v*v) ((vv*)v) (70) 



(and the right-hand-side simplifies even further). We want to show that 

V = (v t v)(w t ) (71) 

but (|70|) in fact only shows that 

V = (v*v)(vv*) + Vi; Viv = (72) 

However, if u is a different eigenvector of A in the same family, that is, if 

Au = fiu (73) 

K[u\ = r m u (74) 

where we assume without loss of generality that fj, ^ A, then the orthogonality @ of u and 
v implies that 

(vv') ((vv*)u) = (v*v) ((vv*)u) (75) 
(since each side is 0). Since the right- hand- side is Vu, we must also have 

Viu = (76) 

so that Vi = when acting on any eigenvector of A (in the same family as v). But this 
means that A and ^4 + Vi have the same decomposition and are therefore equal, which forces 

Vi = 0, and establishes (|7TD . ■ 
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Theorem 2: Let K be the operator constructed via ( \2l\ ) from a given 3x3 octonionic 
Hermitian matrix A, and let y, z G O 3 . Suppose K[y] = ry and K[z] = rz with reK. Then 

(yy*) {{yy ] )z) = {(yy j )z) 

Proof: As in the proof of the Theorem 1, y, z G (T m ) 3 by (j36|), and the components of 
yy^ are in T. We can therefore write 

(di h h\ 
yy ] =\h da h (77) 
\ t 2 h d 3 j 

with d n el and t n G O. Introducing the components y n , z n G T m of y, z, direct computation 
using alternativity shows that 

h(t 3 y 2 ) = 2/3 1 2/i 1 2 1 2/2 1 2 = dihy 2 (78) 
and cyclic permutations, where we have used the fact that 

di = \h\ 2 (79) 
h = ym (80) 

and so forth. But (]46| ) implies that 

t2(hq) = (t 2 t 3 )q+[t 2 ,t 3 ,q}=pq (81) 
for some p G O which is independent of q G T m , and we have therefore shown that 

p = dih (82) 

In particular, we must have 

t 2 (t 3 z 2 ) = d 1 t[z 2 (83) 
and cyclic permutations. A similar argument establishes 

h{hq) = d 2 t 2 q (84) 

for any q G T m , together with its cyclic permutations. 

We are now ready to compute both sides of fl62|) explicitly. We have first of all 



/ dxzi + t 3 z 2 + t 2 z 3 
(yy r )z = t 3Zl + (kz 2 + hz 3 | (85) 
\t 2 £i + hzz + d 3 z 3 

Multiply on the left by yy\ and consider for simplicity only the last component, which is 

t 2 {d 1 z 1 + t 3 z 2 + t^z 3 ) + h(t 3 'z 1 + d 2 z 2 + tiz 3 ) + d 3 (t 2 zx + T x z 2 + ^3^3) 
= (d 1 + d 3 )t 2 z 1 +h(t 3 z 1 ) +t 2 (t 3 z 2 ) + {d 2 + d 3 %z 2 + {\t 1 \ 2 + \t 2 \ 2 + dl)z 3 (86) 
= (d 1 + d 2 + d 3 )(t 2 z 1 + t 1 z 2 + d 3 z 3 ) 
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where we have used flT9|), (|80|), (|3|), and (|84D in the last equality. Repeating this for the 
remaining components, and noticing that 

d 1 + d 2 + d 3 = tr (yyt) = y \ y ( 87 ) 

completes the proof. ■ 
Theorem 1 clearly follows from Theorem 2. But one can also argue that the reverse is 
true, thus avoiding the detailed computation given above. The argument goes as follows: yy^ 
is itself a Hermitian matrix with components in T, and, by the argument given in Section [7], 
it determines the same family structure as the given matrix A, and can therefore be used 
instead of A in applying Theorem 1. 
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